A Solution-Based Adaptive Redistribution Method for Unstructured Meshes

We propose an unstructured mesh redistribution method without using skewed elements for steady-state problems. The regions around solution features are indicated by a sensor function. The medial axes of the strong feature regions are calculated so that elements can be clustered around the most important solution features efficiently. Two approaches, a discrete surface-based approach using a Delaunay triangulation method and a mathematical-representation approach using least square fitting, are shown to calculate the medial axes. Remeshing of an initial volume mesh is performed around the medial axes using an advancing front method and/or an advancing layer method. Two examples are shown to present how our approach works.